Classical Trajectory Diagnosis of Finger-Like Pattern in the Correlated Electron 
Momentum Distribution for Helium Double Ionization 
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With a semiclassical quasistatic model we identify the distinct roles of nuclear Coulomb attraction, 
final state electron repulsion and electron-field interaction in forming the finger-like (or V-shaped) 
pattern in the correlated electron momentum distribution for Helium double ionization [Phys. Rev. 
Lett. 99, 263002; ibid, 263003 (2007)]. The underlying microscopic trajectory configurations respon- 
sible for asymmetric electron energy sharing after electron-electron collision have been uncovered 
and corresponding sub-cycle dynamics are analyzed. The correlation pattern is found to be sensitive 
to the transverse momentum of correlated electrons. 
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Nonsequential double ionization (NSDI) of atoms sub- 
ject to ultrashort intense laser pulses attracts constant 
interests because it is a prototype model for the study of 
three-body Coulomb problem intervened by the highly 
unperturbed interaction of the electrons with the strong 
laser field. Until recently it has been consensus that 
rescattering is the dominant mechanism for NSDI [1| . In 
this three-step mechanism, the first electron is freed by 
a quasi-static tunneling ionization, and is driven back to 
its parent ion and imparts part of energy to dislodge a 
second electron. 

The electron recollision picture as a cornerstone of 
the rescattering mechanism inspires the further investiga- 
tions that achieve insight into the microscopic dynamics 
of the ionization process on the timescale of subfemtosec- 
ond. The advent of experimental techniques, represented 
by the sophisticated Cold Target Recoil Ion Momen- 
tum Spectroscopy (COLTRIMS), combined with high- 
repetition-rate lasers, has to a large extent facilitated 
this type of study. For example, the observed double 
hump structure in the ion momenta Q and the electron 
momenta correlation Q parallel to the field gave solid 
evidence of the time delay introduced by the rescattering 
process and the emission time of both electrons close to 
the zero crossing of the oscillating field. 

Despite the great success of this picture, a comprehen- 
sive understanding of the microscopic dynamics in this 
recollision process is far from being complete. Indeed, 
new high resolution and high statistics COLTRIMS ex- 
periments on double ionization of helium are performed 
independently by two groups and a striking finger-like 
(or V-shaped) structure is observed 0, [B| in the corre- 
lated electron momenta parallel to the laser polarization, 
in qualitative accordance with the prediction of S-matrix 
approach 0, 0, [1] and quantum mechanical calculation 
[g] . These high resolution and high statistics COLTRIMS 
experiments of double ionization (DI) provide benchmark 
data for comprehensive theoretical treatments. 

In this letter, by exploiting an ab initio 3D semiclassi- 



cal model, we have reproduced essentially all the exper- 
imental characteristics, and identified the distinct roles 
of nuclear Coulomb attraction, final state electron re- 
pulsion and electron-field interaction in the formation of 
the finger-like structure. Furthermore, classical trajec- 
tory (CT) analysis facilitates to unveil the sub-cycle mi- 
croscopic dynamics behind the finger-like structure. 

Compared with other approximate approaches exten- 
sively employed in strong field double ionization, e.g. 
the one dimensional quantum model [lOj, strong field S- 
matrix calculation [|| 0, H, [H| and simplified classical 
methods [l2| , our semiclassical model has the advantage 
that all the effects determining the DI ionization pro- 
cess, such as the quantum tunneling, the effective inter- 
actions between particles and the laser field as well as 
the Coulomb focusing effect, can be fully included, while 
keeping the computational capacity still accessible. The 
model has achieved great success in explaining various 
DI phenomena including the excessive DI yield, the 
recoil momentum distribution of doubly ionized ions, mo- 
mentum correlation between two emitted electrons, and 
the energy spectra and angular distribution of photoelec- 
trons. 

For simplicity, we just briefly present the theoretical 
methodology here. We consider a helium atom interact- 
ing with an infrared laser pulse. When the laser field 
is strong enough, one electron is released at the outer 
edge of the suppressed Coulomb potential through quan- 
tum tunneling with a rate given by the ADK formula 
[Til ]. The tunneled electron has a Gaussian-like distri- 
bution on transverse velocity and zero longitudinal ve- 
locity [13]. For the bound electron, the initial position 
and momentum are depicted by single-electron micro- 
canonical distribution (SMD) [15|]. The subsequent evo- 
lution of the two electrons with the above initial con- 
ditions are governed by Newton's equations of motion: 
^p 1 = e{t) — sjn iYne + Vee)- Here the index i denotes the 
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FIG. 1: (Color online), (a) Distribution of correlated electron 
momenta along the laser polarization for Helium DI irradiated 
by 800nm, 4.5 x 10 14 W/cm 2 laser pulses. The black box indi- 
cates the classical limit of the tunneled electron momentum. 
The most important experimental feature, i.e., the finger-like 
structure beyond the limit of 2WUp is well reproduced by our 
model calculation. The model calculations under various cir- 
cumstances yield very different momentum distribution pat- 
terns(see text for details): (b) the laser field is removed and 
the tunneled electrons are replaced by a beam of projectile 
electrons; (c) the electron-electron Coulombic interaction is 
replaced with a Yukawa potential; (d) the nuclear Coulomb 
potential is softened. 



are Coulomb interaction between nucleus and electrons 
and between two electrons, respectively. 

The above Newtonian equations are solved by employ- 
ing the standard 4-5th Runge-Kutta algorithm and DI 
events are identified by energy criterion. In our calcula- 
tions, more than 10 7 weighted (i.e., by the tunneling rate 
) classical two-electron trajectories are traced until one 
electron moves to such a position that > 200a. u.. This 
results in more than 10 4 DI events for statistics. 

The resulting electron momentum distribution, calcu- 
lated with this semiclassical model for the same parame- 
ters as in the experiment [H, is shown in Fig. [TJa). The 
calculation reproduces many key features observed in the 
experiment, including the emission of the two electrons 
primarily into the same hemisphere, the small circular 
accumulation around the zero momentum surrounded 
by four elliptical hard-to-reach regime, and more impor- 
tantly, the finger-like structure beyond the limit of 2^JTT p . 
Here U p = £g/(4aj 2 ), the ponderomotive energy, refers to 
the cycle averaged quiver energy of a free electron in an 
oscillating electric field. 

The off-diagonal and beyond-limit properties of the 
finger-like structure is striking and contradicts to the tra- 
ditional scenario on DI, in which the electron momentum 
at the time of emission is assumed to be small and the 
postcollision electron-electron interaction supposed to be 
weak. Therefore, the parallel momentum 2 of each 



electron results exclusively from the acceleration in the 
optical field: fc[' 2 = smujti on with ionization time 

^ionpiij]- Within this scenario 2-^/T/p should be the max- 
imal momentum and the momentum distribution favors 
accumulation in the diagonal zone because the electrons 
are emitted nearly simultaneously. 

We now proceed to explore the physical effects that 
give rise to this peculiar finger-like structure. In the con- 
text of strong field double ionization, there are essen- 
tially three major effects that may play significant role 
in the double electron emission dynamics: electron-laser 
field interaction that occurs throughout the DI process, 
electron-nuclear Coulomb interaction in the post-collision 
duration and the inter-electron Coulomb repulsion which 
becomes significant when both electrons get close. Be- 
low we investigate all three interactions and clarify their 
distinct roles in the formation of the finger-like structure. 

The first step is to check the role of the external laser 
field and an additional calculation is thus performed, in 
which the laser field is intentionally removed and the 
tunneled electrons are replaced by a beam of projectile 
ones with incident energy of 3.17U P , corresponding to the 
maximal kinetic energy of the tunneled electrons upon 
recollision. The result is shown in Fig. [lib). Two signif- 
icant differences from the complete model calculation in 
Fig. HJa) are found: (i) the finger-like structure beyond 
the limit of completely disappears; (ii) the two 

emitted electrons tend to distribute in the second and 
fourth quadrants of the parallel momentum plane, in- 
dicating that the incident electron transfers much of its 
momentum to the bound one while itself is back-scattered 
into the opposite direction. The comparison between Fig. 
[IJa) and (b) shows the most important role of the laser 
field in turning the two back-to-back emitted electrons 
into the same direction and accordingly the finger-like 
structure. 

The next step concerns the question if this finger-like 
structure is a fingerprint of a strong inter-electron cor- 
relation among the ionizing electrons. We have per- 
formed another calculation in which the final-state elec- 
tron Coulomb repulsion has been deliberately neglected 
by replacing the electron Coulombic interaction V ee = 
with Yukawa repulsion potential of the form V cc ~ 

exp[— Ar;,]/rb, where n, = J | r i t 2 | 2 + b 2 , X = 5.0 and 
b = 0.2. The result of this calculation [Fig. [TJc)] shows 
that the prominent finger-like structure is to a large ex- 
tent reduced and thus provides a clear evidence that the 
final-state electron correlation plays a significant role. 

Last but not the least procedure is to justify the role 
of the electron-nuclear interaction, which is commonly 
believed to be the main reason for the recoil collision 
in field-free (e, 2e) process. This interaction was also 
suggested to be the very ingredient for the field-assisted 
recoil collision in the context of intense field DI of atomic 
helium Q. Accordingly, an additional calculation, in 
which we soften the nuclear Coulomb attraction by em- 
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of interest. 

It has been recognized that the characteristics of the DI 
trajector y ca n be well represented by the recollision and 
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FIG. 2: (Color online), (a-b) Two trajectory configurations 
responsible for the finger like structure. DI yield versus laser 
phase at recollision (c) and at DI moment(d). Here, the statis- 
tics is only collected for the DI trajectories that contribute to 
the prominent finger-like structure, i.e., the regimes of 1.5 
a.u. < |ft{'| < 2.0 a.u. and 2.5 a.u. < \k]j\ < 3.0 a.u., where 
i, j — l,2,i 7^ j. The dashed curves represent the laser field 
for guiding eyes. 



ploying V£ e = — 2/y |rj| + a 2 , where a is chosen as 1.0 
to match the ground state energy of He + , is performed. 

Physically, the shielding of nuclear potential would 
to a great extent diminish the Coulomb focusing effect 
that have significant effects upon both electrons. Clearly, 
a Coulombic potential would attract the tunneled elec- 
tron more dramatically when it moves near the atomic 
core. Such strong attraction may unambiguously bring 
the tunneled electron to share more kinetic energy with 
the bound one. For the bound electron, after achiev- 
ing considerable transferred momentum upon collision, 
it may elastically backscatter from the Coulombic core 
on its way out of the atom. This double scattering pro- 
cess is coined as recoil collision [It} and was routinely 
found in traditional electron impact ionization experi- 
ments, especially when the projectile electron possesses 
the energy of only a few times of the binding energy of 
the inner one. The result shown in Fig. OJd) indicates 
that Coulomb focusing effect is decisive for the produc- 
tion of the electrons with high energy, and thus for the 
finger-like structure. 

The dynamics behind the finger-like structure observed 
in the experiments was supposed to be related with the 
recoil collision in the presence of the external laser field 
This conjecture has been justified by a 2-body 3D 
quantum mechanical simulation of reduced dimensional- 
ity [lH . However the dynamical details of such collision 
is hardly explored in the quantum mechanical treatment. 
In the remaining part of this letter, we proceed to unveil 
the collision dynamics of the DI electrons, especially of 
those that contribute to the finger-like structure beyond 
the limit 2^JTT V . With the CT approach, this becomes 
possible by back-analyzing the history of the DI events 



DI time 12j. We thus provide such an information for 



the trajectories that contribute to the finger-like struc- 
ture in Fig. (2{c) and (d). It is found that, in Fig. HUc), 
the electron pairs contributing to the finger-like structure 
tend to encounter right at zero field. Within rescattering 
picture, this can be understood as that these trajectories 
include the most energetic collisions for which the tun- 
neled electrons are released at the laser phase of about 
17°, and return around the zero crossing of the electric 
field at 270° with maximal energy of 3.17U p [MLIl- 

Upon recollision, the bound electron may be directly 
freed, a process termed as collision ionization (CI), or be 
excited and subsequently ionized by the next field maxi- 
mum, known as collision-excitation ionization (CEI) [20l ]. 
The smaller and larger peaks in Fig. [2]Jd) correspond to 
these two mechanisms, respectively. The smaller ones 
around zero field represents the electron pairs emitted 
after a very short thermalization process (~ attosecond) 
(2l| . While the larger peaks correspond to CEI events 
that are ionized after a few optical cycles delay leading 
to a 0.37T phase difference from the collision phase peaks 
[22I ]. The above calculation indicates that both CI and 
CEI types of DI events contribute to the finger-like struc- 
ture. The time delay in CEI could lead to an off-diagonal 
momentum distribution apart from fcf = k\, but, never- 
theless, cannot account for the excess momentum that 
spills over the limit of 2^JTJ^. 

To further unveil the microscopic mechanism underly- 
ing this unusual pattern, we collect all trajectories that 
constitute the finger-like structure and trace back their 
dynamical evolution at collision and post-collision. With 
the CT diagnosis we sketch two types of trajectory config- 
urations that are responsible for the asymmetric electron 
energy sharing after electron-electron collision and leads 
to the finger- like structure beyond the limit of 2y / [7p (see 
Fig. [21a) and (b), correspond to type-I and II config- 
urations, respectively). While the tunneled electron is 
driven back to the nucleus by the laser field, the field 
strength reduces to zero and the collision is essentially a 
field free three-body system under the pure Coulomb po- 
tential. For type-I trajectories, as shown in Fig. [2)Ja), 
the electron-electron collision near nucleus could lead 
to following consequences: the second electron acquires 
considerable momentum from the returned electron and 
emits in the forward direction, while the returned elec- 
tron is slowed down. Under the influence of the nuclear 
attraction, the latter is transferred to a hyperbolic or- 
bit around the nucleus with a large scattering angle. In 
this way, the returned electron reverses its direction in a 
time scale of attoseconds after the collision. Notice that 
meanwhile the laser field changes its direction. As the 
returned electron has nonzero residual momentum par- 
allel to the instantaneous field direction and is further 
accelerated by the field, its final longitudinal momentum 
is expected to be above the limit 2 A 
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FIG. 3: (Color online). Correlated parallel momentum distri- 
butions with additional conditions on the relative perpendic- 
ular momentum between two electrons, i.e., (a) a.u. ^ ki 2 
sC 0.2 a.u., (b) 0.4 a.u. < k^ 2 sC 0.6 a.u., (c) 1.2 a.u. k^ 2 
^ 1.4 a.u.. (d) The overall relative perpendicular momentum 
distribution. 

electron, after collision its initial momentum is opposite 
to the instantaneous field direction, one expects that its 
final longitudinal momentum is below the limit 
In Fig. [2jb) where a type-II trajectory is schematically 
shown, the situation is similar except that the roles of 
the two electrons have exchanged: Under assistance of 
nuclear Coulomb attraction, the second electron acquires 
a nonzero momentum parallel to the instantaneous field 
and therefore emits with a longitudinal momentum larger 
than 2y / 17p. Although the returned electron is slowed 
down by electron-electron repulsion, it still has a resid- 
ual momentum opposite to the instantaneous field, re- 
sulting in a final longitudinal momentum below 1^JTT p . 
Our statistics reveals that the DI events in the finger-like 
structure consists of 30% type-I configuration and 70% 
type-II configuration. 



The above analysis reveals that the electron-electron 
collision assisted by the nuclear attraction is crucial 
for the emergence of finger-like structure. In certain 
cases that the collision is strong, the finger-like structure 
should be more prominent. To test this idea, we impose 
an additional confinement on our statistics to observe the 
variation of the correlated momentum patterns with re- 
spect to the relative perpendicular momentum between 
two electrons. The results are presented in Fig. [3l 

It is shown in Fig. EKd) that the relative perpendic- 
ular momentum k^ 2 for all DI events distributes over 
an interval of [0,1.4] and exhibits a notable accumula- 
tion around 0.5 a.u.. Moreover, the finger-like pattern is 
more prominent for the case of small relative perpendic- 
ular momentum (see Fig. [3ja)). As we increase the value 
of ki 2 the two finger patterns start to merge (Fig. [3fb)) 
and finally totally disappear (Fig. [3£c)). The above re- 
sult is in agreement with the experimental observations 
[B|,[23|], and provide s a g ood explanation for the ID quan- 
tum calculation in [lfj. Due to dimensional restriction, 
the quantum calculation made there greatly magnifies the 
electron-electron collision effects. As a result, a butterfly- 
like structure similar to Fig. (3Ja) emerges (see Fig. 1 of 

In summary, with a semiclassical quasistatic model we 
have made the CT diagnosis on the finger-like structure 
observed in recent experiments, and unveiled the micro- 
scopic mechanism behind the striking pattern. Our re- 
sults suggest that the finger-like structure is generated by 
the interplay of backscattering of the returning electron 
and the electron-electron Coulomb repulsion. 
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